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Abstract 

We study effective relaxation time of viscous hydrodynamics of strongly coupled non- 
conformal gauge theory plasma using gauge theory/string theory correspondence. We 
compute leading corrections to the conformal plasma relaxation time from the relevant 
deformations due to dim-2 and dim-3 operators. We discuss in details the relaxation 
time r e ff of M = 2* plasma. For a certain choice of masses this theory undergoes 
a phase transition with divergent specific heat cy ~ |1 — T c /T\~ l l 2 . Although the 
bulk viscosity remains finite all the way to the critical temperature, we find that r e ff 
diverges near the critical point as r e ff ~ |1 — T c /T| -1 / 2 . 
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1 Introduction 

Understanding the properties of strongly coupled quark-gluon plasma (sQGP) pro- 
duced at RHIC [1] is an active research area. Of particular interest is the development 
of viscous hydrodynamic codes describing the evolution of the plasma ball at early 
stages of heavy ion collisions [2-4]. These hydrodynamic simulations require input of 
the phenomenological parameters, such as the shear r\ and the bulk ( viscosities, as 
well as higher order transport coefficients necessary to maintain causality of hydro- 
dynamic description [5-8] (and typically, numerical stability of the codes). Currently, 
reliable computation of transport coefficients of sQGP is not available. As a result, one 
resorts to analysis of exactly soluble models of strongly coupled gauge theory plasma. 
Transport properties of a large class of strongly coupled gauge theory plasmas can be 
studied in a framework of gauge theory/string theory correspondence [9,10]. Here, one 
discovers universal relations for the (first order) transport coefficients [11-14] 



where s is the entropy density and c s is the speed of sound in the plasma. The hope is 
that it is generic relations of the type (1.1) that might provide an input into realistic 
simulations of sQGP. 

In this letter we study a specific second-order transport coefficient — the effective 
relaxation time r e ff — in strongly coupled non-conformal gauge theory plasma. The 
motivation for our analysis is the recent work of Song and Heinz [15] where it was 
suggested that the relaxation time near the phase transition can be much larger due 



n = ±. 

S 47T ' 




(1.1) 
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to "critical slowing down". Thus the sQGP produced in heavy ion collisions at RHIC 
might have strong sensitivity to the initial conditions for the bulk viscous pressure. 

We begin in the next section with introducing the notion of the effective relax- 
ation time. The simplest way to obtain non-conformal viscous 4-dimensional hydro- 
dynamics is to compactify a conformal hydrodynamics in (4 + A;)-dimensions on a flat 
A;-dimensional torus [14]. We discuss effective relaxation time of such compactified 
plasma. In section 3 we compute corrections to the effective relaxation time in the 
vicinity of the renormalization group (RG) fixed point induced by relevant operators 
(the mass terms). In section 4 we present results of the detailed analysis of the effective 
relaxation time in M = 2* plasma. We focus on a specific mass deformation of M = 4 
super symmetric Yang-Mills (SYM) plasma which leads to a phase transition in the 
infrared. We summarize our results in section 5. 

2 r e ff from conformal hydrodynamics 

The most general causal viscous relativistic hydrodynamics has many second order 
transport coefficients: 5 - in the case of conformal hydrodynamics [7], and 13 - in 
the case of non-conformal hydrodynamics [16]. It is impractical to simulate hydrody- 
namic evolution in the full multiparameter phenomenological space of these transport 
coefficients. Thus, one typically limits discussion to a single 'effective' second order 
transport coefficient. Since different groups use different approximations of the general 
viscous hydrodynamics, it is necessary to related the analysis through some physical 
observable. We propose to use a dispersion relation of the linearized sound waves in 
plasma to define a common effective relaxation time T e ff. 

u = ±c s k - iTk 2 ± (c* s T eff - k 3 + 0(k 4 ) , (2.1) 

where c s is the speed of the sound waves (obtained from the equation of state), and T 
is the sound wave attenuation (determined by the shear and the bulk viscosities) 

c2 = dv r= (2_5_ i_C_ 

s d£' ' \3£ + V 2£ + V 

The relaxation time r e // is the relaxation time of Muller- Israel-Stewart hydrodynamics 
[5,6]; it is the relaxation time of relativistic conformal hydrodynamics as defined in 
[7]. In general non-conformal hydrodynamics [16], r e // is a combination of particular 
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second-order transport coefficients : 



TV + f£ r n 



^// = - T ^ : - (2-3) 

In A/" = 4 SYM plasma at infinite 't Hooft coupling and in the planar limit, the 
relaxation time is [7, 8] 2 

r eff T ee r^T = . (2.4) 

Starting from the conformal hydrodynamics in d — (4 + A;)-dimensions we can ob- 
tain non-conformal hydrodynamics by compactification on a flat /c-torus [14]. Since 
compactification can not modify the long-wavelength dispersion relation, the hydrody- 
namic dispersion relation for the sound waves (2.1) must be the same. In other words, 
if {ci d \r^,ri^} are the hydrodynamic coefficients of the <i-dimensional conformal 
plasma, then 

where we used an explicit expression for the attenuation in <i-dimensional conformal 
hydrodynamics [7] 

the CFT equation of state £^ = (d — l)V^ d \ and the universality of the shear viscosity 
to the entropy density ratio in strongly coupled plasmas [11-13]. Given (2.2), it is easy 
to see that the bulk viscosity bound proposed in [14] is saturated in such models: 



ci • (2-7) 

V V3 V 

Notice that the effective relaxation time does not change upon compactification. This 
fact, the bulk viscosity ratio (2.7), as well as similar relations for all the other second- 
order transport coefficients of non-conformal hydrodynamics of [16] can be obtained by 
a direct reduction of rf-dimensional conformal hydrodynamics of [7] on a flat A;-torus. 
For example, it can be shown that the "bulk" rn and the "shear" r n relaxation times 
are the same (see also [16]) 

r n = r w = ^ . (2.8) 



^ee [16] for definition of {t w ,tu}- 

2 Finitc coupling/non-planar corrections to the relaxation time of conformal plasmas are discussed 
in [17,18]. 
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The relaxation time of the d — 5 holographic CFT plasma was computed in [19], and 
of the d — 6 holographic CFT plasma in [20] : 




2 1 1 

1 ^ + ^= cotlT 1 ^ In 5 

y/E V5 2 



(2.9) 

2- — -^n3 
Notice that in both cases d 

— >1. (2.10) 

Eq. (2.10) is our first indication that the effective relaxation time in non-conformal 
strongly coupled plasmas is larger than that of the four-dimensional conformal plasma. 
In the following sections we encounter more examples of this phenomenon. 



3 r e ff in the vicinity of the renormalization group fixed point 

An RG fixed point is a conformal theory. The gauge/gravity correspondence [9] not 
only allows us to study a non-trivial (strongly interactive) fixed points (such as Af = 4 
SYM at large 't Hooft coupling), but also provides tools of exploring the deformations 
of such fixed points by relevant (in the infrared) operators (see [21] for a pedagogical 
discussion). In this section we discuss the relaxation time of the strongly coupled 
holographic plasma in the vicinity of the fixed point perturbed by dimension-2 and 
dimension-3 operators. Controllable fixed points that are realized in the holographic 
framework are those of the supersymmetric gauge theories. In these cases we can 
think of the deformation operators as mass terms for the gauge theory bosons and 
fermions. Thermodynamic and hydrodynamic properties of mass deformed conformal 
gauge theories (in the holographic setting) were extensively studied in [14,22-26]. In 
particular, it is straightforward to extend analysis of [24] and compute the 0(k 3 ) 
contribution to the sound wave dispersion relation (2.1). We can then extract the 
effective relaxation time r e //. It is convenient to parametrize the result in terms of the 
deviation parameter 5, 



away from the fixed point: 



T eff = T;(l + p [d] S + 0(5 2 )y (3.2) 



5 



where r* is the universal relaxation time of the holographic conformal hydrodynamics 
at (infinitely) strong coupling (2.4), and (3^ is the correction induced by the operator 
with dim[C] = p. Explicitly we find, 



Again, in both cases the effective relaxation time is larger than that of the conformal 
plasma. 

4 r e ff in J\f = 2* plasma 

We would like now to study the effective relaxation time in a holographic model with 
strongly broken scale invariance. Our choice is the Af = 2* gauge theory at strong 
coupling. Recall that Af = 2* gauge theory is obtained as a supersymmetric mass 
deformation of Af = 4 SYM theory. The low-energy effective description of this theory 
is exactly soluble [27]; furthermore, for a specific point on its Coulomb branch one can 
construct an explicit holographic dual [28] and verify complete agreement between the 
field-theoretic and the dual gravitational descriptions [29,30]. At a finite temperature, 
one has additional freedom of assigning different masses m h ^ rrif to the bosonic and 
the fermionic components of the Af = 2 hypermultiplet [22]. The thermodynamics of 
this theory was studied extensively in [25], and the first order transport coefficients 
in [26]. Here, we extend the analysis of [26] for the dispersion relation of the sound 
waves in Af = 2* plasma to order 0(k 3 ). Using (2.1), we can extract the corresponding 



In what follows, we focus on a specific Af = 2* mass deformation of Af = 4 SYM, 
namely the one with rrif = 0, m& ^ 0. This gauge theory undergoes an interesting 
phase transition for T < T c ps m b /2.29(9), where it becomes unstable with respect to 
energy density fluctuations [25]. Specifically, precisely at T = T c the speed of sound 
waves squared c 2 s vanishes. A perturbative instability of this type is a defining feature 
of a second order phase transition. Since c 2 s oc (T — T c ) 1//2 , the specific heat cy diverges 
as |1 — T c /T|~ 1//2 , suggesting that such a critical point is in the universality class of 
the mean-field tricritical point 3 . Much like real QCD [32], the Af = 2* plasma has a 

3 An identical critical point was recently found in the strongly coupled cascading gauge theory 
plasma [31]. 




p = 2. 



p = 3 



(3.3) 



T eff- 
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Figure 1: (color online) Effective relaxation time r e ff of M = 2* strongly coupled 
plasma. The vertical red line indicates a phase transition with vanishing speed of 
sound. 

plateau in the reduced energy density ^ which extends almost to the phase transition 
- at T ~ \vn,b the Af = 2* equation of state deviates from the conformal equation of 
state by less than 3% [25]. M = 2* plasma has a rapidly growing bulk viscosity in the 
vicinity of the phase transition: at the transition point [14,26], 

c 

v 

which is close to the peak value - ~ 0.7 extracted from lattice QCD 

Techniques for determining the sound wave dispersion relation in strongly coupled 
M = 2* plasma has been developed in [24,26]. We refer the reader for computational 
details to the original work. There is one subtlety though: in analyzing the differential 
equations describing the graviton wave function at order 0(k 2 ) in the momentum (and 
correspondingly its dispersion relation at order 0(k 3 )), we observe that these equations 
become singular in the limit c 2 s — > 0. The latter singularity is removed if instead of 
directly extracting r e ff we compute c 2 r e ff (along with appropriate rescaling of the 
graviton wave function). Upon such rescaling 4 , all the quantities remain finite in the 



C 

= 6.65(3) , or - 

T=T C S 



= 0.52(9), (4.1) 

T=T C 



J A similar procedure has been employed in [14,26]. 
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vicinity of small c 2 . The results of the computations are presented in Figure 1. Given 
that the speed of sound in M = 2* plasma is always less than the conformal value [25], 
i.e., |, we see that for all temperatures down to T c 

Teff > T* . (4.2) 

Moreover, since the quantity c 2 r e ff remains finite at T = T c while the speed of sound 
vanishes as c 2 oc (T — Tc) 1 / 2 , we conclude that r e // actually diverges as T — > T c 

r e// T c oc \1-TJT\- 1 ' 2 , T->T c + 0. (4.3) 



5 Conclusion 

In this letter we studied effective relaxation time in non-conformal strongly coupled 
plasmas using holographic techniques. We discussed theories obtained by compact- 
ification of higher- dimensional conformal plasmas on flat tori, theories obtained by 
deforming a conformal theory with a relevant operator. Finally, using Af = 2* strongly 
coupled plasma as an example, we explored the effective relaxation time in the regime 
with strongly broken scale invariance. In all cases we observed that effective relaxation 
time is longer (and often much longer) than the relaxation time of the M = 4 SYM 
plasma in the planar limit and at infinite 't Hooft coupling [7,8]. In fact, in the Af = 2* 
model we find that the relaxation time diverges near the phase transition with vanish- 
ing speed of sound 5 . This example suggests that the critical slowing down in sQGP is 
a distinct possibility [15]. Further studies of more realistic models of QCD are neces- 
sary to settle this issue. From the holographic perspective, the relaxation time in the 
cascading gauge theory [34] near the deconfinement phase transition should provide a 
useful estimate. 

In the derivation of the second order viscous hydrodynamics from Boltzmann equa- 
tions one finds that the effective relaxation time is [36] 

Boltzmannrp 3?]T T] . 3 . . 

which is much larger than r*T since the ratio of the shear viscosity to the entropy 
density at weak coupling substantially exceed the KSS viscosity bound [12]. This fact, 



5 The same general features of the effective relaxation time can be extracted from the phenomeno- 
logical models of holographic gauge/gravity correspondence studied in [35]. I would like to thank 
Todd Springer for sharing the results of his analysis. 
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supplemented with the analysis of the relaxation time at strong coupling presented in 
this letter, might lead to the speculation that r* is the universal lower bound on the 
relaxation time. Unfortunately, much like it is the case with the KSS bound [37-39], 
this is not so: even though the finite 't Hooft coupling corrections tend to increase 
the relaxation time [17], the non-planar corrections would generically decrease the 
relaxation time [18]. Moreover, in a class of conformal gauge theory plasmas discussed 
in [40], the microscopic causality of the theory imposes the bound 

r eff > 0.56(1) t*. (5.2) 
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